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We solve the quantum constraint equations of the Lemaitre-Tolman-Bondi model in a semiclassi- 
cal approximation in which an expansion is performed with respect to the Planck length. We recover 
in this way the standard expression for the Hawking temperature as well as its first quantum gravi- 
tational correction. We then interpret this correction in terms of the one-loop trace anomaly of the 
energy-momentum tensor and thereby make contact with earlier work on quantum black holes. 
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In the absence of a full quantum theory of gravity, it is 
of interest to consider models which could serve as a pos- 
sible guide in the construction of such a theory [l] . One 
such model is the Lemaitre-Tolman-Bondi (LTB) model 
describing the dynamics of a spherically-symmetric dust 
cloud [2| ■ It has been already used in a variety of papers 
dealing with canonical quantization in both the Wheeler- 
DeWitt framework and loop quantum cosmology, cf. [3j- 
0] and the references therein. While the full quantization 
of the LTB model has not yet been achieved, it was at 
least possible to get insights into the recovery of Hawking 
radiation and black-hole entropy from it. 

Our present paper is a continuation of this earlier work. 
Our motivation is twofold. First, we want to derive a 
quantum gravitational correction to the Hawking tem- 
perature through a semiclassical expansion scheme for 
the quantum states. Second, we want to present an in- 
terpretation of these correction terms through the "trace 
anomaly" of the matter energy-momentum tensor, mak- 
ing thereby a connection to earlier work |8H13| in the 
context of black holes. 

Let us first present the LTB model. The spherical 
gravitational collapse of a dust cloud, in an asymptot- 
ically flat space-time, having energy density e(r, p), is 
described in comoving coordinates (r, p, 9, </>) by the LTB 
metric, 



and 



(d T R) 2 = ^- + 2E=l- F+2E, (3) 

a 

where F(p) = 2GM(p), with M(p) being the active grav- 
itational mass within a p — constant shell, and 



^-i 



(4) 



The function E(p) is the total energy per unit mass 
within the same shell; the marginally bound models are 
defined by E{p) = 0. The case of collapse is described by 
d T R(r, p) < 0. We set c = 1 throughout. 

The canonical quantization of the LTB model was de- 
veloped in [3| and then applied to quantization in a se- 
ries of papers, see [J, [a, 0- Although no full quanti- 
zation has yet been performed, interesting results have 
been obtained at the semiclassical level; they include the 
recovery of Hawking radiation plus greybody corrections 
from solutions to the Wheeler-DeWitt equation and the 
momentum constraints (that is, the quantum constraint 
equations). Insights into the microscopic interpretation 
of black- hole entropy were also obtained, cf. [14 1 and the 
references therein. 

The semiclassical approximation scheme is also em- 
ployed here. We start with the quantum constraint equa- 
tions 0, 



ds 2 



2 (d p R(p,r)) 2 
~ dT + l + 2E(p) 



dp 2 + R 2 (p,r)dn 2 . (i) 



Inserting this expression into the Einstein field equations 
leads to, for vanishing cosmological constant, 



d F 
87rGe(r, p) - " 



R 2 d p R 



(2) 
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where ^> is a functional of the dust variable r(r) as well as 
the gravitational variables R(r) and T(r), and ^p = \/Gh 
is the Planck length. Here, r is the radial variable in the 
ADM formalism ||; we recall that V = F' = 2GM' . We 
note that ([5]) is elliptic outside the horizon and hyperbolic 
inside the horizon; this can be recognized from ((][]). In 
contrast to [J, Q and other papers, no additional factor 
ordering terms are taken into account here; this will be 
crucial in obtaining our results. 



We now make the ansatz 



*[T(r),i?(r),r(r)]=exp(^ r 



drTS(R,r) 



(7) 



where S(R, r) is a function to be determined recursively 
in the following semiclassical approximation scheme. 
With the special factor ordering chosen in [J, Q, this 
ansatz would lead to an exact solution of (|5|) and © 
with semiclassical form. Here, instead, we shall use this 
ansatz to solve ([5]) in a semiclassical approximation; the 
diffeomorphism constraint © is already solved identi- 
cally with this ansatz. 

Inserting (J7J into ([5]), we arrive at 



as 

-ilp6(0) 



r 



-T 



dS_ 
dR 



r 

2T 



d 2 S 



UJU(0).F 



d 2 S 







(8) 
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Here, "(5(0)" indicates the presence of undefined expres- 
sions arising from lim^-^ 6(r — f), which require the pres- 
ence of a regularization scheme. (With the special factor 
ordering chosen in [4j,[5(], the <5(0)-terms are automatically 
cancelled; this is why they were introduced there.) 

For the semiclassical approximation scheme, we make 
the ansatz 

S = S + IpSi + l P S 2 + ... (9) 

and compare consecutive orders in Ip. (The general 
scheme for such an approximation in quantum geometro- 
dynamics is presented in |l|, [15J].) Since T is dimension- 
less, we see from ([7|) that S has the dimension of a length 
(L), so that the dimension of Si is L^ 1 , the dimension 
of S 2 is L~ 3 , and so on. 

Inserting ([5]) into © and comparing different orders in 
lp, we obtain 



■W): [^ 



T 



t)-- 1 - 



which, not surprisingly, is equivalent to the Hamilton- 
Jacobi equation for the action 



A=^JdrTS(R,r) 



(11) 



The next order yields 

,dS dS 1 



oiiV 



r- 



dr dr 
d 2 S Q 



tf 



-i«J(0)- 



dS dSi 
~dR~dR 
d 2 S 



dr 2 



i6(0)F 



dR 2 



0. (12) 



The following order then involves S 2 , 



0(lp) 



dSi 

dr 



TF 
-M(0) 



dS 1 
~dR 

d 2 s 1 

dr 2 



dS dS 2 ~ 




dr dr 

„dS 
1 u 


dS 2 


+ 2 dR 


dR 



W(0) 



d 2 s 1 

dR 2 



0. (13) 



The solution to the Hamilton-Jacobi equation (ITOl) is 



5 = ± ar + / di? 



J" 



+ constant , (14) 



where a = lAv/1 + 2E. With the special factor ordering 
chosen in [J, Q , this would already be the exact solution 
to the Wheeler-DeWitt equation (this, again, was part 
of the motivation to choose that factor ordering in the 
first place). Here, however, the solution (fT4")l occurs at 
the highest order and must be used in the next-order 
equations. 

The variable a above is seen to be related to the (di- 
mensionless) energy E. This is similar to what happens 
for the tunneling mechanism \\S\ Il6j ] . The correspond- 
ing ansatz there involves Lot in place of ar, where t is the 
Schwarzschild time and uj is identified with the conserved 
quantity (in this case, the energy) corresponding to the 
time-like Killing vector. We also note from (|14l) since a 
is dimensionless and So has dimension of L, r has also 
the dimension of L. This will be used later on. 

Observe that the solution (|14l) also holds in the pres- 
ence of a cosmological constant A, with J- then given by 
F = 1 - F/R - AR 2 /3 instead of © , see for the treat- 
ment of positive A and [6j for the treatment of negative A. 
In that case, however, the integral can only be evaluated 
in terms of elliptic functions. Furthermore, the constant 
in (|14p only enters the unknown total normalization for 
^f and will therefore be skipped in what follows. 

Inserting now (|14p into (|12p . we arrive at an equation 
for Si, 



dSi 1 — dSi r . . F 2 - a 2 F 

aT—± + rVl-a 2 F—± + WO— = . 

dr OR R 2 2Fv>l-a 2 F 

(15) 



We solve this equation by the special ansatz 
Si = -Ct-iUi{R) , 



(16) 



where C is a new variable with the dimension L~ 2 be- 
cause Si has dimension L^ 1 and r has dimension L; it 
will play a crucial role below. Since the only length scale 
in this model is GM, we can set for later purpose 



C = 



ai 



(GMf 



(17) 



where a\ is a dimensionless constant. 

Inserting (fl6|) into (fTB"|) . we get a differential equation 
for Ui(R), 



dUi 
AR 



iaC 



VI - a 2 T 



6(0) 



F 



2-a 2 F 



TR 2 2F{1 - a 2 F) 



(18) 



We remark that in the marginal limit, a —¥ 1, this equa- 
tion reads 



dUi _ IR 
~dR- lC ' ! ~= 



6(0) 



F 
F ' 2T(R - F) \ + R 



Integrating (|18|) . we find apart from an irrelevant con- 
stant the desired expression for S± , 



Si{R,t) = ~Ct + aC 
a 2 F 



y/R([l-a 2 }R + a 2 F) 
I -a 2 



(1 - a 2 ) 3 / 2 



In (j{l-a 2 )R+ ^{l-a 2 )R + a 2 F\ 



-i^ (21n(i? - F) - MR + a 2 [F- R\) - lni?) (19) 

We also give the special result for the marginal limit, 
a -¥ 1: 

Si(R,t) = -Ct+ - -i-i^- (2ln(R - F) - lni?) . 

3^/f 2r 

(20) 
In the next order, O (lp), we have to insert the solution 
(HHJ), apart from fI3]). into ([13]). This yields a rather 
complicated equation for S 2 - For example, in the special 
case a = 1, it is of the form 




W(0) fi + i 71 



VT^F 2r R(R-F) 



±2 



R-F dR 



= . 



-3/2 



1(5(0) (g - F) 2 - 2 E 2 
2r R 2 (R-F) 2 



(21) 



In analogy to (|16[) , one could try to solve it with an ansatz 
of the form 



S 2 (R,t) = -Dt-iU 2 (R) , 



(22) 



where D = «2/(GAf) 4 has dimension L~ A and involves 
another dimensionless constant a 2 . We shall not, how- 
ever, follow this here and restrict our attention only to 
the order lp. 

Collecting the solutions up to O (lp) , we can write 

S = So + IpSt = (±a ~ I p C)t ± 



where 



G(R) 



F 
IpaCG(R) - Up5(0)H(R) . 

y/R([l-a 2 ]R + a 2 F) 



(23) 



1-a 2 



a 2 F 



ln(v/(l-a 2 )i? 



(1 - a 2 ) 3 / 2 

+ ^/(l-a 2 )R + a 2 F) , 



(24) 



and 



H(R) = — (2 \n(R - F) - ln{R + a 2 [F-R})~ In R) 



In analogy to earlier papers, cf. [7} , we define positive and 
negative energy states according to the sign in front of the 
dust proper time variable r, with the case of the minus 
sign corresponding to positive energy. Inserting (|23[) into 
the general ansatz (|7|), the positive-energy solution reads 



= cx P[aj /''" r 



"^■^^ 



F 



-IpCr + IpaCG - ilpS(0)H 
while the negative-energy solution is given by 

r-R 



(26) 



<r 



exp (2? ' dr r 

;2 r<~ 1 72 



CLT 



dR 



VI - a 2 F 



F 



IpCr + IpaCG - ilp5(0)H 



(27) 



In order to calculate the Hawking radiation, we shall eval- 
uate the overlap between the "outgoing dust state with 
negative energy" ^^ (where the index "e" refers to "ex- 
panding" cloud) and the "ingoing dust state with posi- 
tive energy" *!>+ (where the index c refers to "collapsing" 
cloud). Since the interpretation of these states is made 
with respect to an observer in the asymptotic regime us- 
ing the Killingtime T, we have to substitute the dust 
time r by T p}. For the outgoing case, we have the 
relation 



T 



dR 



VI - Fa 2 



F 



while in the ingoing case we have 



T = ar - dR 



Vl - Fa 2 



(28) 



(29) 



For the concrete calculation we shall write the full states 
as a product of single-shell states where the radial vari- 
able r is assumed to consist of discrete points separated 
by a distance a. (The continuum limit is obtained for 
a — > 0.) As in @, the Bogolyubov coefficient j3 is cal- 
culated for each shell separately. In the discrete case, 
we replace T by the dimensionless variable 2w and indi- 
cate the dependence on uj by an index. (The factor 2 is 
motivated by the fact that T = 2GM'.) We omit the 
shell index and write the corresponding wave functions 
as ijju (T, R) . We then define j3 to read 



& 



(30) 



where gun is the _R_R-component of the DeWitt metric, 
as it can be read off ((5|) where the inverse of the DeWitt 
metric is lp 4 times the prefactor of the term S 2 ^/5R 2 . 
We thus have gnu = F" 1 ; performing then the required 
coordinate transformation from the variables (r, R) to 
(T, R) gives the result Jg^R = (aF)^ 1 which has to be 



(25) used in the calculation of the Bogolyubov coefficient 



Inserting now ^/ c * and <&+, into (j30|) . we get 



/ dR (glF) -1 exp 



2iwer 
'p 



[T- 



R 



-2iw(7- /" " dR ^ a2-F + 2wa<5(0)iJ ) (31) 



We note that this expression is independent of G. 

In the following, we shall employ a "DeWitt regular- 
ization" and set 6(0) — 0. Whether this can consistently 
be done at the most fundamental level is, however, not 
clear at this stage; here, it is merely used as a formal 
recipe. Recalling that (aJ 7 ) -1 = R/a(R — F), we then 
get 



To evaluate this integral, we use the formula |18| 



which is, in particular, applicable to the case p = and 
< Re v < 1. (We insert a small positive value for Rez/, 
which we let go to zero after the integration.) Using, 
moreover, 

r(-m)T(m) 



i sinh 



u smn ttu 



(with real u), we get 



with 



|& 



2ttF 2 1 
a 2 y e 27T v - 1 



(36) 
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exp 



2iuaT\ 



I 2 
'p 



OC 

F 



exp 



2iu;er 



l + li 



c 



di? 


i? 


i?-F 


Vi 


-a 2 J - 



J" 



(32) 



As in 
able 



we introduce the dimensionless integration vari- 



y 



IujoF 

~1 2 

'p 



i + il- 



di) 



Substituting aus by GAe, 1 where Ae is the energy of a 
shell, and introducing the physical frequency Q = Ae/h, 
we arrive at the final result 



Ift 



2-kGM 



1 



fla 2 (l 



l 2 ^ 

F a 



) «P (*&) 



- 1 



(38) 



and get 



with the quantum-gravity corrected Hawking tempera- 
ture 



exp 



2Fa~ 1 exp 
AiluijF 



2iujaT\ 

) Jo 



I' 2 
tp 



/ 2 
(p 



l + ll- 

a 



' (1 + s) 3 

ds ±— — x 

s 2 +2s 



ds 



(* + l) 2 

s 2 + 2s 



V(l + s) 2 -a 2 (s 2 + 2s) 



(33) 



Up to higher orders of the Planck length squared in the 
exponent, this is so far an exact expression. As in [5J, 
we now assume that the s-integral from zero to infinity 
is dominated by its contribution near s = 0, that is, near 
the horizon; this is also the assumption in the standard 
derivation of the Hawking effect [17| ■ Using therefore in 
P3|) the approximation 



(s + 1) 2 
s 2 + 2s 

2^ V 
we get 

j3 u sa Fa' 1 exp 



v(i 



2- a 



2 -a 2 (s 2 + 2s) 
s + 0{s 2 ) 



(34) 



2kjoT\ 



J 2 
'p 



ds s 



: (i+^ 



exp 



2iusaF 
'p 



1 + / 



C 



(35) 



te^H = 



8^GM(1 + I 2 f) 



(39) 



A number of remarks are in order. First, the meaning 



of the prefactor in the expression for 



(which de- 



pends mildly on fl) is unclear. It is certainly connected 
with the greybody factors, but without a clear-cut nor- 
malization of the quantum states, its interpretation re- 
mains incomplete. Secondly, different from the earlier 
papers, we have calculated the overlap of the quantum 
states in (|30[) for coinciding frequencies us only. The rea- 
son is that an (approximate) thermal spectrum only oc- 
curs in that case. Unlike the highest order Z P , taking 
here also into account two frequencies us and us' corre- 
sponding to two different shells, the integration over us' 
would spoil the thermality. The results (|38p and (|3"9"|) thus 
remain valid only as far as the interaction between dif- 
ferent shells is subdominant. Thirdly, when taking the 
next order in the Planck- mass expansion into account, 
we expect that the term 1 + IpC/a in the denominator 
of (|39[) is augmented by a term proportional to IpD/a, 
where D occurs in (l22l. 



1 Recall that u> is the discretized version of T/2 = GM' , so 
corresponds to GAM = GAe. 



The form of the temperature given in (|39|) for a = 1 
was obtained earlier for the case of the Schwarzschild 
black hole in pjj]. It was calculated there by using the 
quantum tunneling method beyond the semiclassical ap- 
proximation. 

We emphasize that in the previous papers 0, 0] no 
such quantum gravitational correction to the Hawking 
temperature has been found, since calculations have led 
to an exact solution with semiclassical form. 

Substituting now (jTTJ) in (|19p we obtain the expression 
for 5*1 as, 



Si(R,t) 



a i 



(GMf 



Oi\ 



(GMf 
a 2 F 



^R([l-a 2 ]R + a 2 F) 



(1 _ Q a2)3/2 ^ (V(l - a 2 )R + V(l - a 2 )R + ^) 

-i-JfT (21n(i? - F) - ln(R + a 2 [F - R}) - lni?)(40) 

We discuss in the following a method developed in [ll| 
to find the value of the dimensionless constant a\. Con- 
sider for that purpose a constant scale transformation of 
the coefficients of the metric ([1]), given by (cf. also [19j |) 

<W = kg^ . (41) 

Under this transformation we have from ([TJ , 

R = k^R . (42) 

For the Einstein equations to remain invariant under this 
scale transformation, F should according to ((3]) trans- 
form as 



F = k?F , 
and r should transform as 

f = k^T . 

Therefore, (gj| yields 

T = T , 
and since F = 2GM, GM transforms as 



(GM) =k*(GM) 



(43) 



(44) 



(45) 



(46) 



For the Wheeler-DeWitt equation ([5]) to be invariant un- 
der this scale transformation, we must have in addition 
the following transformations: 



f = k-^T , 
where (l4"2"j) and (|4"5")) have been used. 



(47) 
(48) 



Now from the expression for \1/, Eq. (O, we have 
HS^(T,R,t) _ T(r)dS(R,T) 



6X»(r) 



2G dX^{r) 



-*(r,i?,r), (49) 



where /j, = 0, 1 with X a = t and X 1 = R. Because 
we have in (|49p a functional derivative on the left-hand 
side and an ordinary derivative multiplied by T(r) on the 
right-hand side, S is not the action. In fact, as we have 
seen, it is the quantity (fTTj) which has the correct physical 
dimension mass times length and which is equal to the 
action of this model. The first-order action is therefore 
given by 



A 1 = ^JdrTS 1 {R,r) 



(50) 



We have seen that the relevant part for the recovery 
of Hawking radiation was the first, purely r-dependcnt, 
term in (|40)) . Because this part does not contain the grav- 
itational variables, it can for our purpose be considered 
as the matter (dust) action which we shall call A™ . We 
thus have 



A? 



;2 
2G 



dr r 



OtlT 

(GMf 



(51) 



Hence, under the transformations (|44l) . (|4T>|) . and (|48|). 
Af 1 transforms as 



Af 



72 

2G 



drT- 



U\T 



(GM) 2 



k^A'l 1 ~ (1 - 5k)A? , 

where we have assumed that oi\ and r do not scale, and 
k ~ 1 + Sk. Hence, 



SAT = AT - AT = -ATSk 



(53) 



Now using the definition of the energy-momentum ten- 
sor, 



d 4 x, 



2i)A 1 ! 1 



,sk 



dr r- 



(GMf 



2ujcr 



hOL\T 

(GMf 



(54) 



where in the last step, we pass to the discretized version 

by replacing T by the dimensionless variable 2w as before. 

Therefore, by a simple interposition of (|54|) . we obtain, 



O'l 



(GM) 2 
2hujT<j 



-gTH 



(55) 



Next, our task is to find the value of r. For the con- 
tracting cloud case, r is given by the relation (f29|). which 
can be written in a rearranged form: 



r^ + UdR^l 

a a I J- 



(56) 



Now in the calculation of 



A. 



(see equations (|35|) and 



(|36|) b which eventually gives the flux for the model, the 
first term in (|56j) is inconsequential, since it occurs as a 
phase factor (~ e loT ) which yields unity on taking the 
modulus. On the other hand, as we have seen, the terms 
which contribute to the flux come from the integration 
term. Therefore, for the present calculation, considering 
only the last term of ([55]) and using the fact that T = 

1 — C and then substituting s =>/■§ — 1, we obtain, 



t = 2F As 



^^V / (s + l) 2 -a 2 (s 2 + 2 S ) . (57) 



2s 



If this integral is supposed to run from to oo (since R 
runs from F to oo), it is certainly divergent. In order to 
make contact with the expressions in earlier papers, we 
make the following heuristic considerations, which can 
be viewed as a regularization prescription. Restricting 
attention to the relevant regime near s — 0, use of the 
approximation (1341) in the above yields, 



2F 



*4 

2s 



f 



2 _ ° 



= Fins + 0(s). (58) 



Interpreting s as a complex variable and recalling In s 
In \s\ + iargs, we can define a "Euclidean time" te by 



TE 



-7TF 



-2ir(GM) 



(59) 



Finally, substituting this value for te in the Eu- 
clideanized version of (|55|) . we obtain our cherished ex- 
pression, 



OL\ 



GM 
Anhiua 



d x E y 



(60) 



obtain such a correction. Instead, we considered all the 
terms in the expansion for S ^. This led to several 
equations corresponding to different orders in I p. In this 
paper, only terms upto 0(lp) were considered. The equa- 
tions were solved by a special ansatz. After getting the 
solutions for the states upto 0(lp), the "De Witt regu- 
larization" was employed. This regularization enforced 
5(0) — 0. The explicit calculation of the Bogolyubov co- 
efficient near the horizon led to the emission spectrum. 
The corrected Hawking temperature was then automat- 
ically identified. It contained an unknown variable "C" . 
Dimensional arguments then helped us to fix "C" , apart 
from a dimensionless constant. 

The last part of the paper was actually devoted to 
fix the dimensionless constant appearing in "C" . It was 
done by a constant scale transformation of the metric 
coefficients ([1]). A detailed analysis showed that it was 
related to the one loop trace anomaly of the energy - 
momentum tensor for the dust (matter). 

It must be emphasized that a similar result was ob- 
tained by Hawking 19] for an eternal black hole space- 
time by taking into account the one loop correction to the 
partition function due to the fluctuations of the scalar 
fields on the black hole space-time. Exactly the same 
result was also derived later on by different methods 
[g, 0, Ilh - fl3j | . Here our analysis was done in the spirit of 
the quantum tunneling method employing the WKB ap- 
proximation [HH13I ] . Indeed the special ansatz (TTB1 used 
here closely resembles the Hamilton - Jacobi splitting of 
the one particle action S(t,r) = ut + S(r). Such a sim- 
ilarity of our result with Hawking's finding |19J ] may be 
due to the equivalence of the path integral with the WKB 
ansatz upto 0(l p ), a result that has been established ear- 
lier in quite general terms [20( • That this connection also 
holds in the black hole context is a new observation. 



This shows that a\ is related to the "one-loop trace 
anomaly" of the energy-momentum tensor. A similar 
result was obtained earlier in [3, |3, [llH13l | for the eternal 
black hole case. 

To conclude, we mention that a quantum gravitational 
correction to the Hawking temperature from the LTB 
model was established through a semiclassical approxi- 
mation scheme employed in |l0|. Here, no special factor 
ordering [4|, Q was chosen, which was a crucial step to 
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